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ABSTRACT 

We measure the recent star formation history (SFH) across M31 using optical images taken with 
the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit 
the color-magnitude diagrams in ^9000 regions that are ^100 pc x 100 pc in projected size, covering 
a 0.5 square degree area (^380 kpc^, deprojected) in the NE quadrant of M31. We show that the 
SEHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide 
fluctuations in the SEH back to ^500 Myr, most notably in M3Us 10-kpc star-forming ring. We find 
that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ^500 Myr. 

This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. 

We also find that the ring’s position is constant throughout this time, and is stationary at the level 
of 1 km s“^, although there is evidence for broadening of the ring due to diffusion of stars into the 
disk. Based on existing models of M31’s ring features, the lack of evolution in the ring’s position 
makes a collisional ring origin highly unlikely. Besides the well-known 10-kpc ring, we observe two 
other ring-like features. There is an outer ring structure at 15 kpc with concentrated star formation 
starting ^80 Myr ago. The inner ring structure at 5 kpc has a much lower star formation rate (SER) 
and therefore lower contrast against the underlying stellar disk. It was most clearly defined ^200 Myr 
ago, but is much more diffuse today. We find that the global SER has been fairly constant over the last 
^500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SER over the 
past 100 Myr. During the last ^500 Myr, ^60% of all SE occurs in the 10-kpc ring. Einally, we find 
that in the past 100 Myr, the average SER over the PHAT survey area is 0.28 ±0.03 Mq yr“^ with an 
average deprojected intensity of 7.3 x 10“^ Mq yr“^ kpc“^, which yields a total SER of ^0.7 Mq yr“^ 
when extrapolated to the entire area of M31’s disk. This SER is consistent with measurements from 
broadband estimates. 

Subject headings: galaxies: evolution - galaxies: individual (M31) - galaxies: star formation - galaxies: 
stellar content - galaxies: structure 
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1. INTRODUCTION 

A galaxy’s star formation history (SFH) encodes much 
of the physics controlling its evolution. It tells us about 
the evolution of the star formation rate (SFR) through¬ 
out the galaxy, the evolution of the mass and metallic- 
ity distributions, and the movement of stars within the 
galaxy. In addition to the global evolution, focusing on 
the recent SFH (<1 Gyr) reveals the relationships be¬ 
tween stars and the gas and dust from which they form 
and can be used to constrain models of star formation 
(SF) propagation and/or dissolution. For this type of 
study to be possible, however, we first need a spatially- 
resolved view of the past SFH with sufficient resolution 
to probe the relevant physical scales. 

Broad SFH constraints can be derived by looking at 
the properties of the galaxy population across cosmic 
time. However, examining the integrated properties of 
distant galaxies provides limited information on how in¬ 
dividual galaxies form and evolve. While such integrated 
light studies benefit from large sample sizes, the final re¬ 
sults are limited to conclusions about the SFHs of gen¬ 
eral galaxy types (e.g., based on bins of mass, luminosity, 
or color), and cannot say anything definitive about the 
physics that controls the evolution of individual galaxies. 
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To appropriately examine the evolution of individual 
galaxies, it is necessary to study well-resolved nearby 
galaxies for which there exists large amounts of ancillary 
data. With such data, one can, for example, analyze the 
relationship between star SF and gas in the spatially- 
resolved Kennicutt-Schmidt law (e.g., Kennicutt et al. 
2007; Bigiel et al. 2008), understand the evolution of a 
galaxy’s gas reservoir (e.g., Leroy et al. 2008; Schruba 
et al. 2010; Bigiel et al. 2011; Leroy et al. 2013), and 
calibrate SFR indicators (e.g., Calzetti et al. 2007; Li 
et al. 2013), among many others. However, these studies 
have historically been restricted to using only the current 
SFR where ‘current’ is the average over some timescale 
characteristic of a given SFR indicator. These studies, 
therefore, cannot probe the evolution of these relation¬ 
ships with time or on small physical scales where the SFR 
indicators break down (e.g., Leroy et al. 2012). 

For a more detailed analysis of the recent SFH, resolved 
stellar populations are the gold standard. Using individ¬ 
ual stars, we can examine the evolution of a galaxy ar- 
chaeologically by analyzing the color-magnitude diagram 
(CMD) as a function of position within the galaxy. Em¬ 
bedded within the CMD is the history of SF and metal- 
licity evolution of the galaxy. Although recovering this 
information is not completely assumption-free (we must 
make choices about the initial mass function (IMF), stel¬ 
lar models, constancy of SFR within time bins, etc), it 
is the only way to make a time-resolved measurement of 
the SFR. It also has the ability to recover the SFR on 
much finer physical scales. 

CMD fitting has most often been used to probe low 
mass galaxies (e.g., Gallart et al. 1999; Harris & Zaritsky 
2004; Cole et al. 2007; Harris & Zaritsky 2009; Monelli 
et al. 2010; Weisz et al. 2011; Monachesi et al. 2012; Weisz 
et al. 2014) because they are the most numerous type of 
galaxy in the Local Volume, which is one of the few places 
where galaxies can be sufficiently well resolved. The tech¬ 
nique has been used to examine a few individual larger 
galaxies (e.g., Wyder et al. 1998; Hernandez et al. 2000; 
Bertelli & Nasi 2001; Williams 2002, 2003b; Brown et al. 
2006, 2007, 2008; Williams et al. 2009, 2010; Gogarten 
et al. 2010; Bernard et al. 2012, 2015), but these stud¬ 
ies have been limited to either small fields spread across 
the disk and/or halo or have low spatial resolution such 
that it is difficult to pick out detailed features present 
in the galaxy. This technique has never been used to 
contiguously and uniformly recover the recent SFH of an 
L^ galaxy with high resolution. 

In this paper, we present the first finely spatially- 
resolved recent SFH of a significant part of the L^ galaxy, 
Andromeda (M31). M31 has been the target of many 
photometric (e.g.. Brown et al. 2006; Barmby et al. 2006; 
Gordon et al. 2006; Dalcanton et al. 2012; Bernard et al. 
2012; Ford et al. 2013; Sick et al. 2014) and spectro¬ 
scopic (e.g., Ibata et al. 2004; Guhathakurta et al. 2006; 
Kalirai et al. 2006; Koch et al. 2008; Gilbert et al. 2009; 
Dorman et al. 2012; Gilbert et al. 2014) studies due to 
its proximity and similarity to the Milky Way. M31 is 
the ideal place to examine processes in L^ galaxies; it 
is close enough to be well resolved into stars with the 
Hubble Space Telescope (HST) but does not not face the 
same obstacles as studies in the Milky Way, which are 
plagued by uncertainties due to line-of-sight reddening 
and challenging distance measurements. 


We generate maps of the recent (<0.5 Gyr) SFH in 
M31 using resolved stars from recent observations of M31 
taken as part of the Panchromatic Hubble Andromeda 
Treasury (PHAT; Dalcanton et al. 2012). While other 
studies have examined the SFH in M31 using resolved 
stars (Williams 2003b; Brown et al. 2006, 2007, 2008; 
Davidge et al. 2012; Bernard et al. 2012), none have been 
as finely resolved as the work we present here. 

With the resulting spatially- and temporally-resolved 
recent SFHs, we can see where stars form within the 
galaxy and how that SF evolves across the galaxy, 
whether it’s a single star-forming event or propagation 
across the disk. The maps also provide clues about the 
evolution of spatial structure on a variety of different 
scales; while we know a great deal about small-scale SF 
within molecular clouds and large-scale SF within the 
galactic environment, the maps we derive bridge these 
two scales. Recent SFHs also enable the analysis of fluc¬ 
tuations in the recent SFR. This is especially significant 
for SF relations, such as the Kennicutt-Schmidt relation 
(Schmidt 1959; Kennicutt 1989), which often assumes 
a constant SFR over the timescale of the tracer used. 
While this paper deals only with the SFHs themselves, it 
is the first in a series of papers on the SF, dust, and ISM 
contents of M31 on small spatial scales (Lewis et ah, in 
prep). 

This paper is organized as follows: We describe the 
data used in Section 2. In Section 3, we explain the 
method by which we recover the SFHs in each region. 
We present the resulting SFH maps in Section 4 and 
discuss features of the maps in Section 5. We summarize 
the results in Section 6. 

2. PHAT DATA 

We derive the spatially resolved SFHs using photome¬ 
try from the PHAT survey. PHAT surveyed the north¬ 
east quadrant of M31 in six filters, from the near-UV to 
the near-IR, measuring the properties of ^117 million 
stars. Full details of the survey can be found in Dal¬ 
canton et al. (2012) and the photometry is described in 
Williams et al. (2014). Figure 1 shows a 24 /am image 
(Gordon et al. 2006) of M31 with the PHAT footprint 
overlaid. In this paper, we examine the SFH inside the 
solid red region; we have excluded the region closest to 
the bulge (black dashed line) where crowding errors are 
large and the depth of the CMD is shallow, making reli¬ 
able CMD fitting difficult. 

2.1. Photometry and Creation of Single-Brick Catalogs 

We use optical photometry (F475W and F814W filters) 
from the .gst catalogs, which were compiled following 
procedures developed by Dolphin (2000) as described in 
Dalcanton et al. (2012). The stars in this catalog have 
S/N > 4 in both filters and pass stringent goodness-of-fit 
cuts. These cuts leave stars with the highest quality pho¬ 
tometry but with higher incompleteness in more crowded 
regions. The incompleteness is worse in the inner galaxy 
(inside ^3 kpc), which we exclude from our analysis, and 
also affects the centers of stellar clusters. This latter lim¬ 
itation is not a problem for this paper because clusters 
contain only a few percent of the recent SF (Johnson et 
ah, in prep.). Moreover, we are interested in the SFHs 
of the field stars, and picking out stars in the centers of 
dense stellar clusters is not necessary. 
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Figure 1. The PH AT NIR footprint overlaid on a 24 /rm image 
(Gordon et al. 2006) of M31. The thick red line denotes the area 
for which we compute SFHs in this paper. We have left out the 
regions closest to the bulge along the major axis (dashed black 
line) because crowding errors are large and CMD depth is shallow, 
which makes fitting a reliable SFH difficult. The inset shows an 
example of our binning scheme, shown here for Brick 18, which is 
outlined in blue. We have included scale bars for the large image 
as well as for the inset. The bricks are labeled with their numbers 
just exterior to the brick outline. We use these brick numbers 
throughout the text. The image is oriented such that north is up 
and east is to the left. 


The survey was split into 23 regions called ‘bricks’; 
each brick is ^1.5 kpc x 3 kpc in projected size. Odd- 
numbered bricks extend from the galactic center to the 
outer disk along the major axis. Even-numbered bricks 
sit adjacent to them at larger radii along the minor axis. 
Each brick is subdivided into 18 ‘fields’, each with a (pro¬ 
jected) size of ^500 pc x 500 pc. In the optical, adjacent 
fields overlap to cover the ACS chip gap. The SEHs pre¬ 
sented in this paper are derived on a brick-by-brick basis. 
To create a single brick catalog from the 18 ‘field’ cata¬ 
logs, we use the smaller IR brick footprint divided into 
18 non-overlapping regions which roughly describe the 
IR field footprints. In each of these fields, we select all 
of the stars in the corresponding optical catalog that fall 
within the IR field boundaries. We fill in the chip gap 
using two adjacent fields, selecting only the stars that 
fall within the desired portion of the chip gap. Three 
fields in Brick 11 were not observed; this area is com¬ 
pletely covered by overlap from Brick 09 so that there is 
continuous coverage over the survey area. 

The result of this process is the creation of a single 
brick catalog of all stars detected in the optical filters 
that fall within the IR footprint, filling in the chip gap 
and eliminating duplication of stars in the overlap re¬ 
gions. We then grid each brick into 450 approximately 
equal sized, non-overlapping regions that are ^100 pc 
X 100 pc in projected size for a total of ^9000 regions 
across the survey area. In Eigure 1, B18 is outlined in 
blue. The inset shows the binning scheme used within 
each brick. 


Table 1 

50% Completeness Limits and Isochrone Shifts 


Brick 

Number 

^F475W 

(mag) 

mF814lV 

(mag) 

•^log Tgff 


02 

27.2 

26.1 

0.019 

0.22 

04 

27.1 

25.9 

0.019 

0.21 

05 

26.4 

25.0 

0.019 

0.18 

06 

27.2 

26.0 

0.019 

0.22 

07 

26.8 

25.4 

0.019 

0.18 

08 

27.2 

26.0 

0.019 

0.22 

09 

27.1 

25.8 

0.019 

0.21 

10 

27.2 

26.1 

0.019 

0.23 

11 

27.1 

25.9 

0.019 

0.21 

12 

27.3 

26.2 

0.019 

0.24 

13 

27.2 

26.1 

0.019 

0.23 

14 

27.3 

26.2 

0.019 

0.24 

15 

27.4 

26.2 

0.019 

0.25 

16 

27.4 

26.4 

0.019 

0.26 

17 

27.5 

26.4 

0.020 

0.26 

18 

27.8 

26.7 

0.020 

0.29 

19 

27.6 

26.7 

0.020 

0.28 

20 

27.8 

26.9 

0.020 

0.30 

21 

27.8 

26.9 

0.020 

0.30 

22 

27.8 

26.9 

0.020 

0.30 

23 

27.8 

26.9 

0.020 

0.30 


Note. — Column 1 contains the brick number. 
Columns 2 and 3 list the 50% completeness limits in 
F475W and F814W, respectively. Columns 4 and 5 
contain the shifts in log Tgff and used when 

computing the systematic uncertainties. 

2.2. Artificial Star Tests 

Even with the resolution of HST, crowding in regions 
of high stellar density can strongly affect the photometry. 
Many faint stars cannot be resolved in the dense field of 
brighter stars. In addition, faint stars, that would oth¬ 
erwise not be detected, are biased brighter by blending 
with neighboring stars. This also affects brighter stars, 
but to a lesser degree. 

To characterize photometric completeness and to ac¬ 
count for the observational errors that result from crowd¬ 
ing, we perform extensive artificial star tests (ASTs). 
Briefly, we insert fake stars into each image and run 
the photometry as normal. We then test for recovery of 
these fake stars and measure the difference between the 
input and recovered magnitude if a star was detected. 
We adopt the magnitude at which 50% of the stars are 
recovered as our limiting magnitude when solving for the 
SEHs. The completeness limits used for each brick are 
given in Table I. We refer the reader to Dalcanton et al. 
(2012) for further details. 

We inserted ^100,000 artificial stars individually into 
each ACS field-of-view. We combined the resulting ASTs 
into brick-wide catalogs in the same way as the photom¬ 
etry. When running the SEHs for a given 100 pc region, 
we select the fake stars from a 5x5 grid of adjacent re¬ 
gions, such that the ASTs come from a 500x500 pc^ 
region centered on the region of interest. Each of these 
larger regions contains the results of ^50,000 ASTs. 

3. DERIVATION OF THE STAR FORMATION HISTORIES 

We derive SEHs using only the optical data from the 
E475W and E8I4W filters. These filters provide the 
deepest CMDs and the greatest leverage for the recent 
SEHs of interest in this paper. A more detailed discus¬ 
sion of our filter choice can be found in Appendix A. 
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3.1. Fitting the Star Formation History 

We derive the SFHs using the CMD fitting code MATCH 
described in Dolphin (2002). The user specifies desired 
ranges in age, metallicity, distance, and extinction. The 
code also requires a choice of IMF and a binary frac¬ 
tion. It then populates CMDs at each combination of 
age and metallicity, convolved with photometric errors 
and completeness as modeled by ASTs. The individual 
synthetic CMDs are linearly combined to form many pos¬ 
sible SFHs. Each synthetic composite CMD is compared 
with the observed CMD via a Poisson maximum likeli¬ 
hood technique. The synthetic CMD that provides the 
best fit to the observed CMD is taken as the model SFH 
that best describes the data. For full implementation 
details, see Dolphin (2002). 

The fit quality is given by the MATCH fit statistic: fit = 
—2 InL, where L is the Poisson maximum likelihood. We 
estimate the na confidence intervals as > fit — fitmin, 
the Icr confidence interval includes all SFHs in a given 
region with fit — fitmin < 1, the 2cr confidence interval 
includes all SFHs with fit — fitmin < 4, etc. 

We use a fixed distance modulus of 24.47 (Mc- 
Connachie et al. 2005), a binary fraction of 0.35 with 
the mass of the secondary drawn from a uniform distri¬ 
bution, and a Kroupa (2001) IMF. We solve the SFH 
in 34 time bins covering a range in log time (in years) 
from 6.6 to 10.15 with a resolution of 0.1 dex except for 
the range of log(time) = 9.9 - 10.15 which we combine 
into one bin. This time binning scheme was chosen to 
provide as much time resolution as possible while mini¬ 
mizing computing time. We found that using a finer time 
binning scheme with a resolution of 0.05 dex increased 
the computing time by at least a factor of two and only 
resulted in differences in the SFHs of ^1%, which is much 
smaller than systematic and random uncertainties. We 
use the Padova (Marigo et al. 2008) isochrones with up¬ 
dated AGB tracks (Girardi et al. 2010). The [M/H] range 
is [-2.3, 0.1] with a resolution of 0.1 dex. Because we are 
limited by the depth of the data, which does not reach 
the ancient main sequence turnoff, we also require that 
[M/H] only increases with time. We limit the oldest time 
bin to have [M/H] between -2.3 and -0.9 and the youngest 
time bin to have [M/H] between -0.4 and 0.1. 

M31 contains significant amounts of dust (e.g., Wal- 
terbos & Schwering 1987; Draine et al. 2014; Dalcanton 
et al. 2015), which, broadly speaking, can be described 
by three components: a mid-plane component due to ex¬ 
tinction internal to M31 that dominates the older, well- 
mixed stellar populations, a foreground component due 
to Milky Way extinction, and a differential component 
that affects the star-forming regions. In addition to the 
SFH, MATCH allows two free parameters to describe the 
dust distribution: a foreground extinction (Ay) and a 
differential extinction (dAy) which describes the spread 
in extinction values for the stars in each region. The 
differential extinction is a step function starting at Ay 
with a width given by the value of dAy. While fore¬ 
ground extinction is expected to be relatively constant 
across the galaxy, differential extinction can vary signifi¬ 
cantly from region to region as they probe very different 
star-forming and stellar density environments. To de¬ 
termine the best fit to the data, we search extinction 
space to find the combination of Ay and dAy that best 


fits the data. However, the distribution of dust is differ¬ 
ent for young stars and old stars (e.g., Zaritsky 1999). 
The step function differential extinction model provides 
a good fit to the main sequence (MS) component, but it 
cannot reproduce the post-MS stellar populations. 

We mitigate the effects of dust on our SFHs by simpli¬ 
fying the fitting process such that we exclude the redder 
portions of the CMD from the fit. Specifically, we have 
adopted the cuts in Simones et al. (2014), excluding all 
stars with F475W-F814W>1.25 and F475W>21 (shaded 
regions of the CMDs in Figures 2 and 3). This prevents 
contamination from the older populations. We therefore 
avoid extinction-related complications by excluding the 
RGB and the red clump, which is often poorly fit with 
a single step function, and which is not relevant when 
calculating the recent SFH. 

We note that age-metallicity degeneracy is an impor¬ 
tant concern in any kind of SFH work. When model¬ 
ing composite CMDs, it primarily affects the RGB (e.g., 
Gallart et al. 2005), which we do not fit in this analysis. 
Instead, the vast majority of stars in the CMD are main 
sequence stars, for which the age-metallicity degeneracy 
is negligible compared to typical photometric uncertain¬ 
ties. In addition, the metallicity gradient of M31 has 
been extensively studied and found to be very shallow 
(e.g., Blair et al. 1982; Zaritsky et al. 1994; Galarza et al. 
1999; Trundle et al. 2002; Kwitter et al. 2012; Sanders 
et al. 2012; Balick et al. 2013; Lee et al. 2013; Pilyugin 
et al. 2014, Gregersen et al. in prep), and the age range 
we are fitting is small. As a result, we are not concerned 
that the age-metallicity degeneracy affects the results in 
this paper. 

We compute the SFH in 450 regions per brick for 21 of 
the 23 bricks in the PHAT survey. To determine the best- 
fit SFH, we solve multiple SFHs with different combina¬ 
tions of Ay and dAy^ where the best-fit SFH is chosen to 
be the one whose combination of Ay and dAy minimizes 
the fit value as given by the maximum likelihood tech¬ 
nique. Consequently, for each of our regions, we must 
compute many possible SFHs. We minimize the total 
number of SFHs that must be run using an optimization 
scheme to limit the size of (Ay, dAy) space that must 
be searched, as discussed in Appendix B. Based on this 
optimization, we set a constraint that Ay + dAy <2.5. 
In each region, we run a grid of SFHs in (Ay, dAy) space 
with a step size of 0.3 over the range of Ay = [0.0, 1.0], 
also requiring that Ay + dAy <2.5. We take the result¬ 
ing SFH with the best fit, determine the two-sigma range 
around that best (Ay, dAy) pair, and then sample the 
grid in that region down to a finer spacing of 0.1 in Ay 
and dAy. Not only does this ensure that we are finding 
the global minimum, but it also allows us to account for 
the uncertainty in extinction in the SFHs by including 
all fits in the result. In addition, the extinction parame¬ 
ters provide us with an additional method to verify our 
results, as we discuss in Section 3.2. 

As an example, in Figure 2 we plot the CMDs and 
SFHs for two of the regions found near the 10-kpc ring. 
We show the location of these regions over-plotted on 
a GALEX FUV image (Gil de Paz et al. 2007). The 
top region, in Brick 17, is located just off of a spur of 
the 10-kpc ring. The region itself shows very little FUV 
emission, and the resulting SFH is sparse with only mod¬ 
erate SFRs at all times. The lower region falls directly 
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Figure 2. Here we show example results for two different regions, one in Brick 15 sitting on a bright star-forming region (blue, lower) and 
the other in BIT in a more quiescent spot (red, upper). We have overlaid the outlines of Bricks 15 and 17 on top of a GALEX FUV image 
(Gil de Paz et al. 2007). The individual region outlines are also shown. For each region, we show the GMD, with the stars that we fit in 
black and those that we don’t in gray and outlined with red dashes, and the SFH to illustrate the differences on these small scales. The 
SFHs have been binned into 25 Myr increments below 100 Myr and 50 Myr increments above 100 Myr, as the native time resolution (0.1 
dex) allows. Both of these regions have fairly recent SF as can be seen by the well-defined MS in each GMD, though the blue region (B15) 
shows ongoing and more intense SF. Surface densities are calculated using deprojected areas and assuming an inclination of 77° (Brinks & 
Burton 1984). 
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Figure 3. The fit of a region in Brick 15. This is the same region as shown in Figure 2. The left panel shows the data, the middle shows 
the best-fit model for the data, and the right shows the significance of the residuals (the observed CMD minus the modeled CMD, weighted 
by the variance) between the model and the data. All data inside the region outlined by the blue dashed line is masked out of the fit. In 
the left two panels, the color bars indicate the number of stars in each bin. In the right panel, the color bar indicates the significance of 
the residuals. In the residuals panel, darker colors generally indicate more stars in the data than in the model, lighter colors indicate more 
stars in the model than in the data. The lack of systematic residuals in the area that we fit is an indication that the model fits the data 
well. 


on an OB association (OB 54; van den Bergh 1964) in 
Brick 15. As expected, there is elevated, on-going SF in 
that region over at least the past 100 Myr. The CMDs 
for each of these regions show a well-defined MS. The 
consequences of dust are very evident in the region but 
are most easily seen in the part of the of the CMD that 
we do not fit, where the red clump is elongated along the 
reddening vector. 

We note that the Padova isochrones do not include 
tracks younger than 4 Myr. In the resulting SFHs, 
we renormalize the SFR in the youngest time bin to 
reach the present day (0 Myr), conserving the total mass 
formed in that time bin. 

In Figure 3, we show the observed CMD, the best-fit 
modeled CMD, and the significance of the residuals (the 
observed CMD minus the modeled CMD with a weight¬ 
ing determined by the variance) for a region in Brick 15. 
We fit all stars that are outside the blue dashed region, 
which are primarily MS stars, with a smattering of short¬ 
lived blue helium-burning stars. The residuals show no 
distinct features, which means the model is a good fit. 

3.2. Extinction 

In this section we discuss how we incorporate IR-based 
dust maps as a prior in determining our dust parame¬ 
ters. After determining the best-fit SFH in each region 
as described in Section 3.1, we conducted additional ver¬ 
ification by examining the map of total dust. Ay -b dAy. 

We found that there were a handful of regions in which 
the best-fit required large amounts of dust, at or very 
close to the limit of 2.5 mag in spite of there being no ev¬ 
idence for SF within the last 100 Myr, based on a lack of 
luminous MS stars and low SFR averaged over the most 
recent 100 Myr. In these regions, we would expect very 
little dust because there are no dust-enshrouded young 
stars. Bad fits result in these cases because there are 
few stars in the MCD fitting region that can be used to 


constrain the dust. 

To examine this discrepancy, we compared our dust 
parameters with the total dust mass in each region, as 
measured by Draine et al. (2014). We correct the low- 
SFR, high-dust regions by constructing a prior on Ay + 
dAy based on the Draine et al. (2014) dust mass maps 
and multiplying the prior by the likelihood calculated 
by MATCH. The details of the prior are described in Ap¬ 
pendix C. After applying the prior, we compute new fits 
for all of the SFHs measured in each region. We use these 
new fits to go back through Ay, dAy parameter space 
to be sure that we properly sampled 2-cr space around 
the new best-fits for each region. As a result, we were 
able to constrain the MATCH dust parameters in the re¬ 
gions of very-low SFR that are not properly anchored 
in the CMD analysis. Ultimately, because application of 
the prior primarily affects the very-low SFR regions, this 
processing did not significantly affect the SFH results of 
this paper. 

3.3. Uncertainties 

There are three significant sources of uncertainties that 
affect the measured SFHs: random, systematic, and 
dust. In this section, we discuss each source of uncer¬ 
tainty in turn. 

First, we consider random uncertainties. The random 
uncertainties are dominated by the number of stars on 
the CMD and are consequently larger for more sparsely 
populated CMDs. Random uncertainties were calculated 
using a hybrid Monte Carlo process (Duane et al. 1987), 
implemented as described in Dolphin (2013). The result 
of this Markov Chain Monte Carlo routine is a sample of 
10,000 SFHs with density proportional to the probabil¬ 
ity density, i.e., the density of samples is highest near the 
maximum likelihood point. Error bars are calculated by 
identifying the boundaries of the highest-density region 
containing 68% of the samples, corresponding to the per- 
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centage of a normal distribution falling between the ±lcr 
bounds. This procedure provides meaningful uncertain¬ 
ties for time bins in which the best-fitting result indicates 
little or no star formation. 

Next, we consider the systematic uncertainties. Sys¬ 
tematic uncertainties reflect deficiencies in the stellar 
models (i.e., uncertainties due to convection, mass loss, 
rotation, etc.; Conroy 2013) such that different groups 
model these parts of stellar evolution differently, which 
leads to discrepant results for the same data, depending 
on the stellar models used (Gallart et al. 2005; Aparicio 
& Hidalgo 2009; Weisz et al. 2011; Dolphin 2012). These 
uncertainties primarily affect older populations that have 
evolved off the MS. The stellar models of the various 
groups generally agree quite well for MS stars that dom¬ 
inate our adopted fitting region. 

Because we have used the same models across the 
whole survey, all regions experience similar systematic 
effects. We have estimated the size of the systematics 
for a number of regions, covering the range of stellar en¬ 
vironments within M31. We computed the systematic 
uncertainties by running 50 Monte Carlo realizations on 
the best-fit SFHs as described in Dolphin (2002). For 
each run, we shifted the model CMD in log(Teff) and 
Afboi by an amount taken from a random draw from a 
Gaussian with sigma listed in Table 1. These shifts are 
designed to mimic differences in isochrone libraries. We 
then measured the resulting SFH. The range that con¬ 
tained 68% of the distributions from all 50 realizations 
is designated as the systematic uncertainties. The rela¬ 
tive size of the systematics varies greatly from region to 
region but is generally less than half the size of the ran¬ 
dom uncertainties in individual regions and increases at 
larger lookback time. There is also more variation in the 
relative size of the uncertainties in the ring features than 
in the outermost regions where both the stellar density 
and the SFR are low. 

Finally, the variable internal dust content introduces 
uncertainties. We select the best-fit SFH by choosing 
the (Ay, dAy) pair that maximizes the likelihood. How¬ 
ever, there are regions where the difference between the 
fit values of the two most likely SFHs is very small (i.e., 
both SFHs are almost equally likely). We also sample to 
a minimum spacing of only 0.1 in Ay and dAy, and thus 
may have determined a slightly different best-fitting SFH 
than if we had sampled Ay, dAy space more finely. To 
account for these variations, we calculate our uncertain¬ 
ties due to the dust distribution by combining all SFHs 
measured in a given region and determining the range 
that contains 68% of the samples. In this combination, 
the SFHs are weighted by their fit values such that the 
best-fit gets full weight and the na fits are weighted by 
g-o.5n (^e.g., SFHs with fit values that are 2cr from the 
best-fit value are weighted by , or e“^). 

A possible additional source of uncertainty is due to 
the choice of binary fraction, which is a free parameter 
in this analysis. We tested the effect of different binary 
fractions in two of our regions and found that the final 
fit is not very sensitive to binary fraction. This is be¬ 
cause the inclusion of binaries in the model results in a 
color separation on the CMD that is washed out by dust. 
The fits and resulting SFHs are consistent with the un¬ 
certainties when choosing a binary fraction anywhere be¬ 
tween about 0.2 and 0.7. This insensitivity of the SFH 


to binary fraction is consistent with more extensive tests 
presented in Monelli et al. (2010). Uncertainties due to 
binary fraction will be much smaller than those due to 
dust. 

We note that we do not include the model systemat¬ 
ics in our reported uncertainties. While there may be 
absolute uncertainties in the global SFR due to model 
uncertainties, the relative region-to-region uncertainties 
are dominated by the random and dust components. 

3.4. Choice of Region Size 

To generate the spatially-resolved SFH of M31, we di¬ 
vide each of the brick-wide catalogs into regions that are 
approximately 100 pc (projected; 25") on a side, assum¬ 
ing a distance of 783 kpc to M31. There were a few 
different considerations for this size. 

Regions of this size are of scientific interest because 
they bridge the gap between existing knowledge of Galac¬ 
tic pc scale SF (e.g.. Bate 2009; Schruba et al. 2010) and 
SF in more distant galaxies on kpc scales (e.g., Leroy 
et al. 2008). The resolution is also fine enough to resolve 
features such as large HII regions and giant molecular 
clouds. 

While a finer grid would also be scientifically interest¬ 
ing, there are a couple of difficulties to consider. The 
main problem is that smaller regions would have insuf¬ 
ficiently populated CMDs, increasing the random uncer¬ 
tainties of the SFHs to unacceptable levels. With our 
adopted ^100 pc bin size, the number of stars within the 
GMD fitting region ranges from 110 to 3900. In 93% of 
the regions, we fit more than 500 stars, and in 80% we fit 
more than 1000 stars. Additionally, the SFHs are com¬ 
putationally expensive to run. For each region, MATCH 
must be run multiple times to determine the (Ay, dAy) 
pair that provides the best fit to the data. Moving to 
50 pc size regions would have resulted in four times as 
many regions, significantly increasing the time needed to 
derive the SFHs. Even at our 100 pc grid size, deriving 
the SFHs and uncertainties for the entire sample required 
more than 500,000 GPU hours using XSEDE resources 
(Towns et al. 2014). 

Our overall technique is similar to that of Simones et al. 
(2014), who measured the SEHs of UV-bright regions 
within Brick 15 of the PH AT survey. Their goal was 
to convert the SEHs into EUV fluxes and compare with 
the observed fluxes in each star-forming region. About 
half of their regions had fewer than 500 stars on the part 
of the GMD they were fitting but they still found rea¬ 
sonable agreement between the modeled and observed 
fluxes. This agreement indicates that the GMD fitting 
routine is robust, even with a modest number of stars in 
the part of the CMD occupied by young stars. 

3.5. Reliability of the SFHs as a Function of Lookback 

Time 

We have excluded the red side of the CMD in our SEH 
recovery process, so consequently, our fits are not sensi¬ 
tive to old stellar populations. The exact age at which 
sensitivity is lost is set by the oldest stars observable on 
the MS, which varies with stellar density and dust ex¬ 
tinction. 

We perform two different tests to examine the sensitiv¬ 
ity of our results as a function of time. Eirst, we create 
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Figure 4. Artificial CMDs over a range of extinction and pho¬ 
tometric depth, assuming a constant SFH, solar metallicity, and 
three different extinction values; one low: [Ay, dAy] = [0.0,0.0], 
one typical of our SFHs: [Ay, dAy] = [0.3,0.8], and one high: [Ay, 
dAy] = [0.6,1.8], shown in the left, middle, and right columns, re¬ 
spectively and labeled in the upper right corner of the top panel of 
plots. The brick number listed in the left panel of plots indicates 
the depth of the corresponding row. Within each panel, individ¬ 
ual stars are color-coded by age. All stars older than 1 Gyr are 
dark red. In the top 3 rows, we have over-plotted 630 Myr, so¬ 
lar metallicity isochrones, one attenuated by Ay, and the second 
attenuated by Ay + dAy. Individual stars within the CMD will 
be reddened to somewhere between the two isochrones. Increased 
extinction reddens some stars off of the MS and into the region we 
do not fit in the SFH recovery process (shaded gray). The bottom 
two rows each have isochrones of 4 different ages over-plotted. The 
ages are listed in the left panel. 


artificial CMDs using MATCH. The CMDs are generated 
with a constant SFH and solar metallicity, while mod¬ 
eling observational uncertainties by using the results of 
the ASTs in each region. The results are shown in Fig¬ 
ure 4, where we plot simulated CMDs. The youngest 
stars are shades of blue, and all stars older than 1 Gyr 
are red. Each of the top three rows shows a CMD at 
a different depth, where the top row is the shallowest 
CMD closest to the bulge (B05), and the third row is 
the deepest in B23. The brick numbers are indicated 
in the left panel. The columns display varying amounts 
of extinction, which is applied to the CMD in the same 
way we apply it to the model CMDs when recovering 
our SFHs. Extinction is labeled in the top panel by [Ay, 
dAy]. The left column is un-reddened, the middle col¬ 
umn shows the effect of the median extinction found in 
each of our regions ([Ay, dAy] = [0.3, 0.8]), and the right 
column shows the upper limit of extinction allowed in 
our SEHs. We also check the sensitivity by over-plotting 


isochrones of a single age (630 Myr) and solar metallicity. 
The isochrone in the first column has not been reddened. 
In the second and third columns of the first three rows, 
we plot two isochrones, one extincted by Ay and one ex- 
tincted by Ay + dAy. Individual stars can be extincted 
to anywhere between the two isochrones. 

The variation in depth is large across the survey. The 
inner regions are very crowded and as a result the pho¬ 
tometric depth is quite shallow. In Brick 05, where there 
is high crowding and moderate levels of extinction, we 
cannot detect many stars that are older than ^500 Myr. 
As we move further away from the center of the galaxy, 
crowding and extinction generally decrease and we can 
probe stars that are 700 - 1000 Myr old (row 3 of Eigure 
4). The exception, of course, is the 10-kpc ring, where 
extinction can be high; the CMDs in the second row of 
Eigure 4 mimic the conditions found in these regions. At 
low or mid-range extinction, we still see many stars on 
the MS that are at least 600 - 700 Myr old. As we in¬ 
crease the extinction, many of these stars are reddened 
into the portion of the CMD that we do not fit, as can 
also be seen from the isochrones. Unreddened 630 Myr 
old stars are easily recovered; however, stars that have 
the highest level of extinction applied to them are en¬ 
tirely reddened into the neglected portion of the CMD. 

In the last two rows, the colored CMD is identical to 
that in the second row (B15 depth). Here, though, we 
plot solar metallicity isochrones of different ages: the 
fourth row has isochorones of 100, 200, 300, and 400 Myr 
while the bottom row shows isochrones of 500, 630, 800, 
and 1000 Myr. Erom these tests, we can see that there 
will be a number of regions that are reliable back to 1 
Gyr, while others that are more extincted may only have 
a handful of stars that are more than ^500 Myr old. 

Ultimately, due to the wide variety of stellar environ¬ 
ments found across the survey, not all bricks can reliably 
cover the same age range. The inner regions are much 
more crowded and have shallower completeness limits 
than the outer regions; as a result, the age limit for re¬ 
gions closest to the bulge is only 400 Myr. Regions that 
fall within the ring and spiral arm features are much more 
extincted than other regions; although the completeness 
limit is deeper, stars may be reddened into the region 
of the CMD that we do not fit. These regions have an 
age limit of no more than 500 Myr. In the outer regions 
where extinction and crowding are less, the completeness 
limit is much deeper and the SEHs are reliable back to 
700 - 1000 Myr. We therefore choose a conservative age 
limit of 400 Myr for consistency across the survey. Eor 
all scientific analysis, we examine only the most recent 
400 Myr. In some of the plots that follow, we include 
results back to 630 Myr for illustrative purposes, with 
the caveat that they are only relevant for the outermost 
regions of the survey. 

4. RESULTS 

We now present the results of the SEH fitting process. 

4.1. Star Formation Rate Maps 

We combine the single best-fit SEHs in each region 
into maps of Esfr as a function of time. As an ensem¬ 
ble, they reveal the recent SEH of M31 in the PH AT 
footprint over the last 630 Myr. This SEH is shown in 
Eigure 5. Each panel displays the SER surface density 
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Figure 5. Map of the SFH of M31 covered by the PHAT survey. Each region is colored according to its SFR surface density in each time 
bin. We have applied a lower cut at a Ssfr of 10“"^ M© yr“^ kpc“^ in order to highlight structure in the image. The time range covered 
is shown in the lower left of each plot. The plot in the top left shows the GALEX FUV image, smoothed to the same physical scale as 
our SFHs for better comparison. A blue dashed curve is over-plotted on each panel to aid the eye in recognizing structural change between 
timebins. The maps are oriented as in Figure 1. Note: We have included the time bins from 400 Myr to 600 Myr for illustrative purposes 
only. As discussed in Section 3.5, the SFHs in these time bins are much less certain in the crowded inner regions and in the dusty ring 
than they are in the outer parts of the disk. 
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Figure 6. SFR surface density as a function of radius and time. 
Here we plot the SFR surface density in all regions within a 40 
degree arc about the major axis (shaded yellow in the inset) per 
time bin. We have binned the radius in 0.5 kpc bins. The SFR 
surface density in each radius bin is the mean of all regions that 
fall within that bin. The thick black line shows the overall mean. 
The dashed line shows the fractional RGB residuals (Seth et ah, in 
prep.) with values given on the right axis. 

within the time range specified in the lower left corner. 
We note that these time bins are not the native resolu¬ 
tion of A(logt) = 0.1; instead, we have binned the results 
within the most recent 100 Myr into ^25 Myr bins. This 
helps to expose the continuous structure, especially in 
the most recent 25 Myr when SF is very low across the 
galaxy. The regions are colored according to Ssfr in 
that region during the given time range. To more clearly 
illuminate the structure, we have placed a lower limit on 
the SFRs visible within this map such that all regions 
with SsFR lower than 10“^ Mq yr“^ kpc~^ are colored 
black. We have also over-plotted a blue dashed line on 
each image to aid the eye in recognizing structural evolu¬ 
tion between time bins. There is little large-scale change 
over the last ^500 Myr. The upper left panel shows 
the GALEX FUV image smoothed to the same physical 
scale as our SFHs, which shows excellent morphological 
agreement with the most recent time bins. This agree¬ 
ment strengthens confidence in our SFH maps given that 
measurements made in ^9000 completely independent 
regions reproduce coherent large-scale structure seen in 
a well-established SFR tracer. The relationship between 
the SFR and the FUV flux will be examined in an up¬ 
coming paper (Simones et al. in prep.). 

The SFHs reveal large-scale, long-lasting, coherent 
structures in the M31 disk. There are three star-forming 
ring-like features; a modest inner ring at ^5 kpc, the 
well-known 10-kpc ring, and an outer, low-intensity ring 
at ^15 kpc that partially merges with the 10-kpc ring 
due to a combination of projection effects and a possible 
warp that is visible in HI (e.g.. Brinks & Shane 1984; 
Chemin et al. 2009; Corbelli et al. 2010). These rings 
have been observed previously in /Infrared Ar¬ 

ray Camera (IRAC) images (Barmby et al. 2006) and 
Spzter/Multiband Imaging Photometer (MIPS) images 
(Gordon et al. 2006), as well as in atomic (Brinks & 
Shane 1984) and molecular gas images (Nieten et al. 
2006) . There is also an observed over-density of red giant 
branch (RGB) stars in the 10-kpc ring (Dalcanton et al. 


2012). Recovery of these coherent, well-known features 
is further confirmation that our method is robust. 

One of the most remarkable features of this SFH is 
that the 10-kpc ring is visible and actively forming stars 
throughout the past ^500 Myr. Although SF occurs at a 
low level in the outer regions at all times, SF in the ring 
feature at 15 kpc is most concentrated starting at 80 Myr 
ago. This is likely due to SF in OB 102, which has had an 
elevated SFR compared to its surroundings over the last 
^100 Myr (Williams 2003a). The inner ring feature at 
^5 kpc is also visible, and though there appears to be SF 
in that ring distinct from the surrounding populations, 
that feature gains definition 200 Myr ago but has largely 
dispersed in the last 25 Myr. 

We further investigate these trends in Figure 6 where 
we plot the average SFR surface density as a function of 
radius in each time bin for a subset of the regions along 
the major axis. We determine the distance to the center 
of each region and then divide the regions into bins of 0.5 
kpc. The three rings at ^5 kpc, ^10 kpc, and ^15 kpc 
are clearly visible as peaks in Hsfr, providing further 
indication of ongoing (if low) SF in the ring features over 
500 Myr. This trend is further supported by evidence 
of RGB star residuals in the ring features (Seth et ah, 
in prep ), as shown by the thick dashed line. There 
is an over density of RGB stars in the 10-kpc ring and 
in the inner ring where our oldest time bin shows the 
highest SF. The plot reveals that not only is the 10-kpc 
ring long-lived, it has also remained mostly stationary 
in galactocentric radius over 500 Myr, a result that has 
implications for the origin of the ring, which we discuss 
further in Section 5.4. 

4.2. Mass Maps 

In Figure 7, we show the evolution of recent mass 
growth in the galaxy from 630 Myr to the present. The 
upper left panel shows the 3.6 jiva image, which is a rough 
estimation of the total mass of a galaxy, smoothed to the 
same spatial resolution as our SFHs. The next 11 pan¬ 
els show the total mass formed over a given time range 
(the same as those in Figure 5). The upper middle panel 
shows the mass that formed in the last 25 Myr while the 
bottom right panel shows all of the mass formed in the 
last 630 Myr. The time range is written in the lower left 
of each plot. In the upper right, we have indicated the 
total mass formed and the average SFR over that time 
range. 

Most of the SF in M31 occurred much earlier than 
the timescale probed by these SFHs; consequently, the 
amount of mass accumulated over the last ^500 Myr has 
been minimal. As can be seen in the upper left image in 
Figure 7, the older stars are distributed fairly uniformly. 
This means that the structure that we see in our inte¬ 
grated mass maps only appears within the last Gyr 
(more recent than the ages of stars that dominate the 
emission at 3.6 jam). Over the time range of our SFHs, 
SF has been confined primarily to the rings. 

We further examine the impact of SF in the last ^500 
Myr by looking at the fraction of mass formed during this 
time. In Figure 8, we plot the fraction of mass formed 
in the last 400 Myr to the total mass in each region. We 
derive the total stellar mass by multiplying the 3.6 jam 
image (Barmby et al. 2006) by a constant mass-to-light 
ratio of 0.6 (Meidt et al. 2014). Only ^0.8% of the total 
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Figure 7. Map of the cumulative mass formed. Each region is colored according to the mass formed in that region, integrated from the 
present day, and scaled by the area of the region. The time range covered is shown in the lower left of each plot. In the upper right of each 
plot, we indicate the SFR surface density and the total mass formed over that time range. We have chosen to present the maps in this way 
in order to better highlight the structure of the most recent mass formed. Most of the mass in M31 was in place by at least 1 Gyr ago. We 
scale the maps to the mass formed in the last 630 Myr in order to highlight structure. The maps are oriented as in Figure 1. Note: We 
have included the time bins from 400 Myr to 600 Myr for illustrative purposes only. As discussed in Section 3.5, the SFHs in these time 
bins are much less certain in the crowded inner regions and in the dusty ring than they are in the outer parts of the disk. 
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Figure 8. Map of the fraction of mass formed in the last 400 
Myr compared to the total mass as inferred from the 3.6 iim image 
(Barmby et al. 2006). The amount of mass formed over the time 
range covered by these SFHs is very small, as expected, with the 
highest fractions coming in the 10 and 15-kpc ring features. The 
map is oriented as in Figure 1. 

stellar mass of this section of the galaxy formed in the last 
400 Myr, which is ^3% of a Hubble time. On a region- 
by-region basis, the mass fraction reaches a maximum of 
^7% and is highest in the ring features, specifically in 
the two outer features, where most of the gas is located 
(see Figure 9). Consequently, though the SFR in the ring 
in the last ^500 Myr is much higher than that in the rest 
of the galaxy, it must still be very small relative to past 
SFRs. 

5. DISCUSSION 

We have presented maps of SF and mass evolution in 
M31 that show rich structure with ongoing SF and evo¬ 
lution. Much of this structure has been observed in maps 
of M31 in other tracers. As an example, in Figure 9, we 
plot our SFR map averaged over the last 100 Myr, next 
to maps of FUV flux (Gil de Paz et al. 2007) and HI 
(Brinks & Shane 1984). All three maps show M31’s ring 
structures. In addition, we note the good agreement be¬ 
tween the 100 Myr SFR map and the map of FUV flux, 
which is sensitive to SF within the last 100 Myr. These 
regions of high SFR and flux are also broadly consistent 
with areas of largest gas content. We now discuss some 
of the features observed in these maps. 

5.1. The Recent PH AT SFH 

The SFHs of the individual regions are extremely di¬ 
verse. Some of the regions, such as those in the 10-kpc 
ring have SF that is ongoing and seemingly long-lived; 
others, such as those that fall in the outer parts of the 
galaxy, have generally quiescent SFHs in the past 500 
Myr. At the native time resolution of 0.1 dex, the un¬ 
certainties on the individual SFHs in each region are sig¬ 
nificant. By combining all the SFHs, we can reduce the 
amplitude of the uncertainties in those time bins and 
obtain a more significant and constrained result for the 
total SFH. 

We derive the total SFH within the PHAT footprint 


Table 2 

Total PHAT SFH 


log(ti) 

log(yr) 

log(tf) 

log(yr) 

SFR 

Moyr-i 

6.60 

6.70 

0-09l°;“ 

6.70 

6.80 


6.80 
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6.90 

7.00 


7.00 

7.10 

n OO+0.02 
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7.80 

7.90 
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7.90 
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8.50 

8.60 

0 39+0-01 


Note. — The total SFH 
summed over all regions. The 
first column lists the start of 
each timebin and the second 
column lists the end of each 
time bin. The third column 
shows the total SFR in each 
time bin, integrated over the 
entire survey. To convert to 
SFR surface density (Ssfr), 
divide the SFR by the total 
area of the survey, 378 kpc^. 

The uncertainties represent the 
smallest range that contains 
68% of the probability distribu¬ 
tion calculated fom the random 
and dust uncertainties. The 
first time bin extends to the 
present day. 

by integrating over all regions in Figure 5. In Figure 
10 we show the total SFR per time bin over the survey 
area. Figure 11 shows the same SFH but with linear time 
bins over 3 different time ranges. The dashed blue line 
indicates the average SFR over the past 100 Myr. The 
uncertainties in each of these figures include the random 
component as well as uncertainties in the Ay, dAy com¬ 
bination of the best-fit SFH of each region. There are 
also systematic uncertainties due to isochrone mismatch 
(see Section 3.3), which we do not include but are ^30% 
in all time bins. The SFRs and uncertainties in each time 
bin are listed in Table 2. 

We look more closely at the total SFH to determine 
which regions contribute in each time bin. While the 
SFR has been generally declining over the last 600 Myr, 
we see a bump at ^50 Myr. Further examination reveals 
that this peak is not a galaxy-wide feature. To show this, 
in Figure 12, we plot the SFH in 5 angular slices in theta 
about the major axis. We see that the slice that falls 
down the major axis (green line) shows a declining SFR 
over all periods of time, i.e., there is no 50 Myr bump. 
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Figure 9. Maps of the average SFR in the last 100 Myr as derived in this paper, GALEX FUV (Gil de Paz et al. 2007, 4.3" resolution), 
and HI (Brinks & Shane 1984, 45" resolution). There is good morphological agreement between the SFR map and the FUV map, which 
is sensitive to the same time scale. It is also clear that the regions of highest SF, primarily in the 10-kpc ring, coincide with the regions of 
highest HI content. Orientation of the images is as in Figure 1. 


The west-most slice (blue line) shows a slight peak at 50 
Myr and if we move to the slice that is out along the mi¬ 
nor axis (purple line), the SFR shows the most prominent 
50 Myr bump. These slices contain bright OB associa¬ 
tions and are dominated by ring regions. Consequently, 
it is likely that the 50 Myr peak in the total SFH is due 
to specific star-forming regions in the rings. 

In Figure 10, we have over-plotted a dashed line show¬ 
ing the average SFR over the past 100 Myr. The av¬ 
erage SFR in the PHAT area over the past 100 Myr is 
0.28T0.03 Mq yr“^, where the uncertainty is calculated 
from the tests of systematic uncertainties. To determine 
the average SFR over the past 100 Myr over the entire 
disk, we adopt a simple scaling argument. We use both 
the FUV and 24-/im maps of M31 and select the D25 
ellipse (Gil de Paz et al. 2007) as an aperture. We then 
measure the total flux within the entire aperture and the 
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Figure 10. The total SFH of M3I within the PHAT footprint 
over the last 400 Myr, combining the individual SFHs of each field. 
The dashed line (blue) shows the average SFR density over the 
most recent 100 Myr. The red error bars are a combination of the 
random uncertainties and the uncertainties in Ay and dAy. The 
time bins are A(logt) = O.I, so the youngest bins cover consider¬ 
ably less linear time than do the older bins. 


flux inside the PHAT footprint (without Bricks 1 or 3). 
In both cases, the fraction of the total flux inside the 
PHAT footprint is ^40%. Therefore, to scale the PHAT 
SFR to the entire D25 aperture, we must multiply by a 
factor of 1/0.4 = 2.5. We multiply the SFR by this scal¬ 
ing factor to yield a total SFR of ^0.7 Mq yr“^. Previ¬ 
ous studies have examined the global SFR using resolved 
stars (^1 Mq yr“^; Williams 2003b), 8 jam (0.4 Mq yr“^; 
Barmby et al. 2006), FUV (0.6-0.7 Mq yr“^; Kang 
et al. 2009), Ha (^0.3 Mq yr“^;Tabat abaci & Berkhui- 
jsen 2010, 0.44 Mq yr“^; Azimlu et al. 2011), and FUV -h 
24 jam (0.25 Mq yr“^; Ford et al. 2013). Our result falls 
within this range, and is most consistent with methods 
that effectively average over longer timescales. 

Although this simple scaling argument neglects possi¬ 
ble variations in the SFR from one side of the disk to 
the other, we do not expect these variations to be large 
enough to dramatically affect the overall SFR. A com¬ 
plete analysis of the SFR measured via different tracers 
is beyond the scope of this paper but will be examined 
more closely in an upcoming paper (Lewis et ah, in prep). 

5.2. Birthrate Parameter 

We further examine the total SFH by computing the 
birthrate parameter, 6, which is defined as the ratio of 
the current SFR to the past averaged SFR: 


We calculated the past-averaged SFR by converting 
the 3.6 /im image to mass, as described in Section 4.2 
and setting the lifetime of the disk to the Hubble time. 

While the birthrate parameter is generally used to clas¬ 
sify a galaxy as a whole as in either a star-bursting or 
quiescent phase, in this case we can examine b for our in¬ 
dividual regions. This gives us the ability to tell whether 
regions of current high SF, such as the 10-kpc ring, are 
currently undergoing unusually high SF relative to the 
rest of the galaxy (in which case b should be large) or 
whether those regions have always had higher SFRs (in 
which case b < 1). In Figure 13, we show a map of the 
birthrate parameter over two timescales: the lifetime of 
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Figure 11. Same as in Figure 10 but showing time bins on a linear scale. The first panel shows the SFH back to 400 Myr. The second 
shows the SFH over the last 100 Myr and the third panel shows the SFH within the last 25 Myr. 


our SFHs (^400 Myr) and 100 Myr. The regions are 
colored according to 6 on a linear scale. We see over 
both time periods that the rings at 10 and 15 kpc have 
higher values of b than the rest of the galaxy. While b is 
greater than unity over most of these ring features, there 
are very few regions where it rises above 2, which indi¬ 
cates that while SF is certainly elevated in those regions 
compared to the overall SFR, it would not be considered 
‘bursty’. The exception in the right panel of Figure 13 is 
the large OB associations in Bricks 15 and 21 where we 
know that SF is occurring at an elevated rate. These low 
values of b are not wholly unexpected since the 3.6 /im 
image is quite smooth (see the upper left panel in Figure 
7) though one can still make out the ring feature. 

Integrated over the entire survey area, we find that the 
birthrate parameter is 0.23 over the last 100 Myr and 0.27 
over the last 400 Myr. This further shows that M31 has 
been forming stars in the last 400 Myr at a much lower 
rate than it was in the past. 

5.3. Comparison with Previous Work 

Here we examine the results presented in this paper in 
relation to earlier studies of the recent SFH of M31. 

Williams (2002) measured the SFH using archival HST 
data in various fields throughout the galaxy, some of 
which overlap with the PH AT survey area. He also recov¬ 
ered a SFH using CMD analysis, with an earlier version 
of the code used in this work. The overlapping fields in 
his study fall primarily along the 10-kpc ring, so they 
probe regions of higher SFR. His ‘INNER’ field, which 
falls in B13, shows a SFH with a steep decline from 1 Gyr 
to ^200 Myr ago, then a short rise to a peak in a time 
bin covering the range 40-80 Myr ago, and a decrease to 
the present day. Two nearby fields, G287 and Gll, show 
a similar morphology in their SFHs, though the peak 
comes at slightly earlier times, closer to 100 Myr. Other 
fields fall along the northeast section of the 10-kpc ring 
and show strong SF at the most recent times. Our SFHs 
are consistent with this. As we showed in Figure 12, the 
average SFR per time bin can vary significantly in differ¬ 
ent regions. If we were to define our regions differently, 
such that they fell only on the regions of strongest SF, 
we would also see a sharper rise to the present day. 


Williams (2003b) measured the SFH of M31 using 
ground-based data (Massey et al. 2006) and found an 
increase in the SFR from ^25 Myr ago to the present. 
This increase was seen primarily in the northeast spiral 
arm. While the work we present here shows a decrease in 
the SFR from ^50 Myr to the present, our last time bin 
in Figure 5 from 25 Myr to the present shows a tighten¬ 
ing of the ring structure and more localized SF along the 
northeast arm, suggesting that the difference in resolu¬ 
tion between this study and that of Williams (2003b) 
could be the cause of the discrepancy. For example, 
blends could be measured as upper MS stars, artificially 
increasing the measurements in the youngest time bins 
from ground-based data. 

Davidge et al. (2012) examined the recent SFH of the 
entire disk of M31 by comparing u' luminosity functions 
with those derived from models assuming various SFHs. 
Their result indicates a factor of 2-3 rise in the SFR 
during the past 10 Myr, which is in broad agreement 
with the results of Williams (2003b). 

While our results do not show a recent SFR increase, 
we will point out that Williams (2003b) and Davidge 
et al. (2012) looked at the SFH over the entire disk of 
M31, while we focus only on ^1/3 of the disk. There is 
evidence that the SFR is elevated in the outer regions 
of the southern and western parts of the disk more so 
than in the eastern parts. Our study does not include 
the southern disk, nor does it reach the outer regions of 
the western edge of the disk. It is possible that these 
regions could contribute to an increase in the last 10 - 
25 Myr that we do not see in this work. 

Williams (2003b) also saw movement of SF across the 
disk, which he interpreted as evidence of propagating 
density waves from the northern to the southern disk. 
We do not see this movement. However, our survey does 
not cover the southern disk and our region sizes are much 
smaller, giving us the ability to more precisely locate SF. 
We find no evidence for propagation in the area covered 
by the PH AT survey. We do, however, agree that SF 
has been confined primarily to the ring structures over 
at least the last 500 Myr. 
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Figure 12. SFH in slices about the major axis. On the left we show the average SFR per area in each time bin for all of M31 divided 
into 5 angular slices. The lines are colored according to angular slice. The same coloring is used in the right panel to indicate the angular 
slice that is represented. All slices cover 30 degrees, though they do not all contain the same number of regions. The black line shows the 
average SFR per time bin for the entire survey area. 



Figure 13. Maps of the birthrate parameter b over 400 Myr and 100 Myr timescales. The regions are colored according to the value of 
the birthrate parameter in that region and an upper limit has been set at 6 = 2. The maps are oriented as in Figure 1. Both maps show 
6 > 1 in the 10 and 15-kpc rings. The OB associations in Bricks 15 and 21 are clearly visible in the right image as regions of elevated b. 
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Figure 14. SFH showing the average SFR in each time bin for all 
regions within the 10-kpc ring (red) and all other regions (black). 
Error bars include random uncertainties primarily due to the num¬ 
ber of stars on the CMD as well as uncertainties as a result of the 
search in Ay and dAy for the best-fit SFH. Comparison with Fig¬ 
ure 10 shows that the total SFH is driven by SF within the ring. 
SF outside of the ring feature was more important before ~500 
Myr ago. It is important to note that due to the high inclination 
of M31, our definition of the 10-kpc ring, shown in the inset, may 
also include some of the 15-kpc ring (if it is indeed a distinct fea¬ 
ture). Nevertheless, this still confirms our result that SF occurs 
primarily in these ring features. 

5.4. The Mystery of the 10-kpe Ring 

5.4.1. The Ring is Long Lived 

The results we present here indicate that the ring of 
SF at a radius of 10 kpc has persisted for at least 500 
Myr. If we assume a rotational velocity of M31 at the 
10-kpc ring of r\j 250 km s ^ (Chemin et al. 2009; Corbelli 
et al. 2010), then the dynamical time of M31 at the ring 
is ^250 Myr. We have defined the dynamical time as 
the time to make one full rotation at a given radius: 
^dyn = 27rrIVr. This means that SF has continued in the 
ring for at least two dynamical times. Williams (2003b) 
found that SF has occurred in the ring for the last 250 
Myr, and Davidge et al. (2012) found SF in the ring for 
at least 100 Myr. Dalcanton et al. (2012) showed that 
there is an over density of stars with ages >1 Gyr in the 
ring as well, further supporting the long-lived nature of 
this feature. 

Looking at just the regions that fall inside the 10-kpc 
ring, we see that over the last 600 Myr, SF in the ring 
drives the overall SFR (see Figure 14). Within the last 
400 Myr, the SFR in all regions that do not he in the 
ring has been mostly constant, while the SFR inside the 
ring has varied substantially. This is also clearly visible 
in Figure 5, where the ring is very prominent due to a 
lack of SF occurring in other parts of the galaxy. We find 
that in the last 400 Myr, ^60% of all SF occurs in the 
10-kpc ring feature. 

The long-lived nature of the ring is also visible in Fig¬ 
ure 6, where SF is elevated relative to the areas just 
outside of the ring at all time. Only in the most recent 
25 Myr do we see a notable decrease in SFR in the 10- 
kpc ring. This plot only includes regions that fall along 
the major axis so as to minimize uncertainties in the 


deprojected distance to each region. As a result, this 
slice does not include many of the higher SFR OB as¬ 
sociations, such as OB 54 and those that fall along the 
northeast portion of the 10-kpc ring. These would likely 
act to increase the SFR over the most recent ^50 Myr. 

Not only is the ring long-lived, but it has also been 
mostly stationary. From Figure 6, we can see that the 
10-kpc ring has remained centered at about the same 
location, moving no more than 0.5 kpc over ^500 Myr. 
This translates into a motion of ^1 km/s. 

5.4.2. Dispersion of the Ring 

Our maps of SF show that the 10-kpc ring is clearly 
broader at older ages than at younger ages. This suggests 
two possibilities: 1) Over time, SF in the ring has grown 
more concentrated, occurring in a much narrower strip, 
likely following the density of gas within the ring features 
or 2) SF has always occurred in the center of the ring and 
we see a broadening of the ring as the stars disperse. The 
second option is more probable, since we know qualita¬ 
tively that the majority of stars form in clusters and that 
those clusters eventually disperse and the individual stars 
diffuse into the surrounding environment, becoming part 
of the larger galactic background (e.g., Harris & Zaritsky 
1999; Bastian et al. 2009). 

The simple kinematical argument is that if stars are 
born with an average random velocity of 10 km s“^ and 
that motion is directed radially outward, then over the 
course of 100 Myr, those stars should move 1 kpc. For the 
results we present here, this would result in significant 
motion over the course of 500 Myr (5 kpc). It would 
also suggest a strong dispersion of the ring features over 
relatively short periods of time, which we do not see. In 
reality, these timescales are much longer. Bastian et al. 
(2011) showed that the evolution of the spatial structure 
of stellar populations in nearby dwarf galaxies occurs on 
many different timescales, setting lower limits of tens to 
hundreds of Myr, longer than would be expected based 
on simple arguments. 

A detailed quantitative discussion of the evolution of 
spatial structure will be the subject of a future paper. 

5.4.3. What is the Origin of the Ring? 

The 10-kpc ring is the most prominent feature of M31, 
visible in wave bands from the UV through the radio 
(e.g.. Brinks & Shane 1984; Barmby et al. 2006; Beaton 
et al. 2007; Block et al. 2006; Gordon et al. 2006; Massey 
et al. 2006; Nieten et al. 2006; Gil de Paz et al. 2007; 
Braun et al. 2009; Ghemin et al. 2009; Fritz et al. 2012). 
Its ubiquity naturally leads to the question of its origin. 
Surprisingly, there have been few studies devoted to this 
question; its genesis remains largely uncertain. 

Those studies that have attempted to answer the ques¬ 
tion have invoked one of two likely formation scenarios 
for a ring: (1) resonance ring due to a central bar, or (2) 
collisional ring due to the passage of a satellite galaxy 
through the disk of M31. 

Bars can produce inner, outer, and/or nuclear rings 
in galaxies (Buta 1986). M31 shows evidence of a bar, 
as suggested by its boxy bulge seen in infrared imaging 
(Beaton et al. 2007), though it is likely not very strong 
(Athanassoula & Beaton 2006). The bar length has been 
estimated to be 4 - 5 kpc (Athanassoula & Beaton 2006; 
Beaton et al. 2007). If the 10-kpc ring is due to a bar. 
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it is unlikely to be either a nuclear ring, which generally 
occurs within the bar (e.g., Buta & Combes 1996) or an 
inner ring, originating at the end of the bar (Schwarz 
1984). Instead, if the ring is indeed due to a rotating 
bar, it must be an outer ring, which occurs near the 
outer Lindblad resonance (OLR, e.g., Schwarz 1981; Buta 
1995). It has also been shown that in barred galaxies 
with rings, the ratio of the outer ring diameter to the 
bar diameter is ^2 (Kormendy 1979; Athanassoula et al. 
2009). By comparing observational data of M31 with 
N-body simulations, Athanassoula & Beaton (2006) find 
that the OLR is at 45 ± 4 arcmin (9 - 10 kpc at the 
distance of M31). The outer ring radius to bar length 
ratio and the location of the OLR in M31 suggest that 
the 10-kpc ring could indeed be due to a bar. 

Rings in galaxies can also be the result of collisions 
with satellite galaxies (e.g., Lynds & Toomre 1976). In 
the case of M31, that satellite galaxy is likely M32. A 
handful of studies have attempted to model the effect of 
M32 crashing through the disk of M31 (Gordon et al. 
2006; Block et al. 2006; Dierickx et al. 2014) and have 
found that such a collision could produce a ring with 
properties similar to that observed, including ring centers 
offset from the center of the galaxy, and the hole seen in 
the 10-kpc ring. These studies suggest that the impact 
event happened 20, 210, and 800 Myr ago, respectively. 

There are inconsistencies with both scenarios. It is 
much more statistically likely that M31 is a barred galaxy 
with resonance rings simply because there are many 
barred galaxies that show this kind of ring structure. 
However, the offset centers of the 5 and 10-kpc rings 
(Block et al. 2006) and the hole in the 10-kpc ring are 
difficult to explain if the rings are all due to a bar. On 
the other hand, the collisional simulations suffer from un¬ 
certainties in the mass and orbit of M32, both of which 
strongly affect the collision and resulting morphological 
features. In addition, if a collision occurred, we would 
expect to see evidence of ring expansion. We note, how¬ 
ever, that it is possible that existing spiral structure and 
bar disturbances could interact strongly with these much 
weaker collisional rings generated by a companion, pre¬ 
venting the expected classical propagation (Struck 2010). 

Our results indicate that the ring has been a dis¬ 
tinct feature, actively forming stars for at least 500 Myr. 
Davidge et al. (2012) found that the last major distur¬ 
bance to the disk occurred at least 500 Myr ago. Bernard 
et al. (2015) and Williams et al. (2015) both find the last 
major event in M31’s history to be 2 - 3 Gyr ago. These 
are longer than the timescales suggested by Gordon et al. 
(2006) and Block et al. (2006) for a collision with M32. 
We cannot rule out a collision that happened 800 Myr 
ago (Dierickx et al. 2014) much less 2-3 Gyr ago. How¬ 
ever, the velocity deviations suggested by collisional ring 
models do not match the observations presented here. 
Block et al. (2006) expect the rings to expand radially 
at ^7 - 10 km s“^ and Dierickx et al. (2014) find radial 
motions of at least 20 km s“^ in their models. The sys¬ 
tems described by Struck (2010) have radial motions of 
at least 10 - 20 km s“^, though he does note that M31 
seems to be an exception to the general class of colliding 
ring galaxies. 

Consequently, based on existing models of the origin of 
M31’s ring features, we can rule out a purely collisional 
origin because of the additional time-resolved data we 


present in this paper. A bar-induced ring is more likely 
to be long-lived and resonance rings from bars are fairly 
common. If the ring is a resonance ring, it would make 
M31 a much more average galaxy than one that had re¬ 
cently suffered a dramatic collision. On the other hand, 
perhaps the morphological features we see in M31 result 
from a combination of the two scenarios. We note that 
no existing studies attempt to explain the existence of 
the outer ring at 15 kpc. 

If we wish to fully understand the origin of the 10-kpc 
ring (or either of the other ring features) in M31, it is 
clear that more work is needed. One necessary step is 
to run simulations that follow the creation of the ring 
and its evolution over 1 Gyr or more. These simulations 
should account for multiple possible formation scenar¬ 
ios: collision with a satellite, resonance due to a bar, 
and a combination of the two. Athanassoula & Beaton 
(2006) provide stellar velocity field predictions that need 
to be tested. Finally, obtaining deep photometry that 
goes below the oldest main sequence turnoff over a large 
area around the ring in combination with a more com¬ 
plex dust model would allow us to determine the SFH in 
M31 back to more than 5 Gyr ago, providing additional 
observational constraints on the lifetime of the ring. 

6. CONCLUSIONS 

We have measured the SFH of ^1/3 of the star-forming 
disk of M31 from HST images as part of the PHAT sur¬ 
vey. We divided the survey area into ^9000 approxi¬ 
mately equal-sized regions and determined the SFH, as 
well as the extinction distribution, independently in each 
region. 

We have found that SF in M31 has been largely con¬ 
fined to three ring features, including the well-known 10- 
kpc ring, over the past ^500 Myr. The 15-kpc ring be¬ 
comes most prominent and structured starting at about 
80 Myr ago, with SF increasing to the present day. In 
the 5-kpc ring, SF reached a peak about 100 Myr ago 
and the ring feature has since begun to dissolve. 

The 10-kpc ring is long-lived and stationary, producing 
stars over at least the past 400 Myr, which is about two 
dynamical times at that radius. Over this period of time, 
SF has been elevated relative to the surrounding regions. 
As a result, the ring drives the overall SFH. The shape 
of the total SFH follows that of the SFH of just regions 
that fall broadly in the 10-kpc ring feature. This feature 
has also shown little or no significant propagation in our 
survey area over the lifetime of these SFHs. 

We have shown that the total mass formed over the 
last 400 Myr is less than 10% of the total mass of the 
galaxy, as measured from 3.6 fim images. M31 is a fairly 
quiescent galaxy today. Aside from a handful of bright 
regions that coincide with OB associations, the total SFR 
at the present day is very low. We have shown with 
maps of the birthrate parameter, 5, that only a handful 
of regions have b > 2 over the last 100 Myr or 400 Myr. In 
most cases, the current average SFR is very low relative 
to the past average. Globally, the galaxy has b = 0.23 
over the last 100 Myr and b = 0.27 over the last 400 Myr, 
further confirming that M31 is a fairly quiescent galaxy 
today. 

Finally, aside from a possible increase in SF ^50 Myr 
ago, the SFH over the last 400 Myr has been relatively 
constant. We have computed the average SFR over the 
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past 100 Myr to be 0.28±0.03 Mq yr“^ within the PHAT 
footprint. Extrapolating to the entire galaxy we find a 
global SFR over the last 100 Myr of 0.7 Mq yr“^, which 
is consistent with the range found in previous studies. 
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APPENDIX 
A. FILTER CHOICE 

We have chosen to use only the optical filters (F475W and F814W) for this analysis. In the PHAT survey, data 
were also taken in two NUV filters (F275W and F336W) and two NIR filters (FllOW and F160W). The optical filters 
provide the greatest leverage for probing the recent SFH because they have by far the deepest CMDs of MS stars 
which allows us to probe both further down the luminosity function and further back in lookback time. 

In principle, addition of the UV filters might potentially allow for greater constraints on the recent SFH; however, 
in practice, the PHAT data do not allow it. Only the brightest main sequence stars have measurements in the two 
UV filters, which significantly reduces the age range over which the SFH can be derived. As such, including the UV 
filters will not improve the analysis presented in this paper. 

The NIR data are more complex. They are significantly shallower than the optical data. The mean depth in F160W 
is an absolute magnitude of ^0, which is 2 magnitudes shallower than in the optical. At this depth, we lose the main 
sequence at less than 200 Myr. In addition, there are not very many stars that fall on the main sequence to begin 
with, further reducing the usefulness of the NIR for analysis of the recent SFH. Finally, the NIR models of young stars 
are much more uncertain than the optical models, as there has been less verification of the models in this regime. 

In Figure 15, we plot the color-magnitude diagrams of a bright, highly star-forming region in B15. The left panel 
shows the optical data and the right panel shows the NIR data. On each panel, we have plotted three sets of isochrones 
at 25, 100, and 400 Myr at two different metallcities, solar and slightly sub-solar ([M/H] = -0.3). M31 has a very flat 
metallicity gradient that is approximately solar, so these values are reasonable for the disk. The stars mark the main 
sequence turnoff points at each age and metallicity. The optical data are deep, and we have leverage back at least 400 
Myr. On the other hand, the IR data are much shallower and do not allow reliable SFH measurements beyond ^200 
Myr. Consequently, we have chosen to use just the optical data in this analysis of the recent SFH. 
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Figure 15. The optical (left) and NIR (right) CMDs of a highly star-forming region in BIS. On each CMD, we have plotted three sets 
of isochrones at 25, 100, and 400 Myr and two metallcities: solar ([M/H] = 0.0, red) and sub-solar ([M/H] = -0.3, blue). We also show a 
reddening vector of 2.5 magnitudes which is the amount of extinction required by the best-fit SFH in this region. In the left panel, we also 
show the optical exclude region described in this paper (dashed line), the 50% completeness limit (dot-dash line). 
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B. DETERMINING THE BEST-EIT SEH BY SEARCHING Ay, dAy SPACE 

The two main free parameters in the fitting process are those that describe the dust model: an extinction applied 
evenly to all stars, Ay^ and a differential extinction, dAy^ that is applied in addition to Ay. The differential extinction 
is applied as a distribution to all stars, rather than uniformly, so that each star receives some additional extinction 
between 0 and dAy. As a result, the total extinction applied to each star falls in the range [Ay, Ay + dAy]. 

Using B15, we sampled Ay, dAy space on the range Ay = [0.0, 1.0] and dAy = [0.0, 6.0] with step sizes of 0.25 
in both parameters. The results for two of these regions are shown in the top panel of Figure 16. We found that 
while some regions had a well-defined best-fit SFH, other regions showed “troughs” of (Ay, dAy) pairs going to 
arbitrarily high dAy values with fits that were within 1-cr of the best-fit, defined as /itbest + 1- The resulting SFHs 
are indistinguishable within the error bars. These “troughs” are an artifact of the exclusion area we used in the SFH 
solution, which removed the red clump and RGB from the fitting process. At higher extinction, stars are artificially 
pushed into the exclude region on the CMD. In order to account for the fewer stars in a given CMD bin, the SFR 
increases. Examination of the regions with these troughs shows that they tend to start at Ay + dAy = 2.5. As a 
result, we limit our total extinction to 2.5 magnitudes. 
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Figure 16. Exploration of Ay^ dAy parameter space in two regions of B15. The top panel show examples of a grid search in Ay^ dAy 
space. One region has a very well defined best fit. In the second region, a “trough” of low fit values occurs. If dAy is increased, the fit 
value does not change, but rather remains very close to the best fit. In the bottom panel, we show the results of the search over (Ay, 
dAy) space for the same two regions with the requirement that Ay + dAy <2.5. Once again, the region on the left remains well-defined 
in Ay, dAy space. In the second region, the constraint on extinction has forced a best fit. The lines show the I-, 2-, and S-cr contours. 



































22 


Lewis et al. 


150 


300 


NmS stars 

450 600 


750 

I 


900 


1050 



SFRxIO- 4 [M© yr- 


Figure 17. Comparison of the MATCH-derived dust parameter {Av+ dAy) with the dust mass surface density as derived by Draine et al. 
(2014) before (left) and after (right) applying the prior. On the left axis, we show the ratio of the MATCH extinction to the dust mass surface 
density. On the right axis we show the ratio of the MATCH extinction to the extinction derived from the dust mass surface density following 
Equation Cl. The relationship between the two parameters is generally quite constant with an expected increase in scatter at the lowest 
SFRs. 


C. APPLYING A DUST PRIOR 


In some of the very low-SFR regions at the survey edges and between the ring features, we found that MATCH assigned 
large extinctions with very poor constraints. Given the low SFR and the small number of bright, young MS stars at 
these locations, there is no physical explanation for large extinction. To investigate this issue, we use the Draine et al. 
(2014) dust mass maps to examine the relationship between the MATCH extinction and the total dust mass in those 
regions. On the left side of Figure 17, we plot the ratio of the MATCH total dust {Ay + dAy) to the dust mass surface 
density as a function of SFR (averaged over the last 100 Myr) and number of MS stars. Each circle represents a 
single region, color-coded by the number of MS stars in that region, where we define MS stars as those with F475W — 
F814W < 1 and F475W < 26. The bricks closer to the bulge have lower completeness, so the value of 26 was chosen to 
accommodate those bricks with 50% completeness just fainter than 26th magnitude. On the right axis, we converted 
the dust mass surface density into extinction according to Equation 7 in Draine et al. (2014): 


Ay = 0.74 


^Md 


10^ Mq kpc 


30 


(Cl) 


We have made this conversion for ease of reader interpretation, though analysis of this result is beyond the scope of 
this paper. Eor the sake of this exercise, only the relative numbers are important. 

If there were a perfect correlation between measured dust mass and the MATCH dust parameters, we would expect 
to see a straight horizontal line in this figure such that the ratio of the normalized MATCH dust to the normalized dust 
mass would be constant across all SERs with some scatter expected at the low SER end. Instead, we see that at very 
low SERs there is wide scatter in this ratio, and, as may be expected, these regions also contain lower numbers of MS 
stars. If there are very few or no stars on the upper MS, it is more difficult to anchor the SEH. In order to fit the 
broadened lower MS, the amount of differential reddening is increased. Eigure 17 shows that the ratio between MATCH 
dust model parameters and the dust mass is approximately constant in regions with high SER and/or many MS stars. 
This implies that we should be able to fit a straight line to a plot of MATCH dust parameters vs. dust mass in these 
regions with the most robust fits. So we use this relation to set a prior on the total extinction in a region to adjust 
the low-SER, high-dust regions to achieve more physical extinction parameters. 

We compare the MATCH best-fit extinction {Ay -b dAy) with dust mass surface density (M^ust) in each region. To 
do this, we fit a line with slope m and scatter a to a plot of MATCH dust {Ay + dAy) vs. M^ust of the n regions with 
high SER and many MS stars. 

The slope and dispersion of that line are found by solving the following equation numerically for m and a such that 
X" = l: 
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is the number of stars 


where Ay + dAy is that of the best-fit SFH for that region and = C/\/]V^tars- stars 

in the nth region and C is calculated such that the resulting values of m and a do not depend on A^stars- fhis case, 
C = 12. We split into two solutions to account for the fact that we have set an upper limit of Ay + dAy = 2.5 in 
our analysis. As we have shown in Appendix B, the best-fit extinction parameters are effectively lower limits in some 
of our regions. As a result, we assume that all regions with Ay + dAy > 2 are lower limits. 

We find that a line with 777 = 4.5 x 10“^ and a = 0.6 results in the fit that best helps us constrain the regions with 
low SFR and high dust. We apply this prior to our results in all regions by recomputing the fit values for each {Ay, 
dAy) pair. The new fit value is given by 

I + dAv) - m X Mdust 
jT-tnew — f'^told T ^2 

We compare the Ay and dAy values corresponding to the new fits in each region with the dust mass. The results are 
shown on the right side of Figure 17. Applying the prior did tighten up the relation a little bit, reducing the number 
of outliers on both tails at the low SFR end. 











